UNCLASSIFIED 


AD  -  4  2  5  3  0  1 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


I 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any- 
patented  invention  that  may  in  any  way  be  related 
thereto. 


AFCRl-6a-39‘ 
AO  425  301 


ADMITTANCE 
OF  A  WAVEGUIDE 
RADIATING  INTO 
STRATIFIED 
PLASMA 

SCIENTIFIC  REPORT  NO.  2 
CONTRACT  N0.AF19|628)-2410 

AUGUST  1963 


At#  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
OFFICE  OF  AEROSPACE  RESEARCH 
LAURENCE  G.  HANSCOM  FIELD 
BEDFORD,  MASSACHUSETTS 

APPLIED  RESEARCH  LABORATORY 


'SYIAAX I  A  ELECTRONIC  SYSTEMS 

( invcrmncnt  Systems  Management 

GENERAL  TELEPHONE  &  ELECTRONICS^, 


AFCRL-63-394 


ADMITTANCE  OF  A  WAVEGUIDE 
RADIATING  INTO  STRATIFIED  PLASMA 

By 

Jan is  Galejs 


AFPLIED  RESEARCH  LABORATORY 
3 YL VANIA  ELECTRONIC  SYSTEMS 
A  Division  of  Sylvania  Electric  Products  Inc. 
40  SYLVAN  ROAD.  WALTHAM  54,  MASSACHUSETTS 


August  1963 


Contract  No.  AF1 9 ( 628) - 2410 

Project  No.  _ 4642 _ 

Task  No.  464202 


SCIENTIFIC  REPORT  NO.  2 


Air  Force  Cambridge  Research  Laboratories 
Office  of  Aerospace  Research 
Laurence  G.  Hanscora  Field 
Bedford,  Massachusetts 


NOTICES 


Requests  for  additional  copies  by  Agencies  of  the 
Department  of  Defense,  their  contractors,  and  other 
Government  agencies  should  be  directed  to  the: 


DEFENSE  DOCUMENTATION  CENTER  (DDC) 
CAMERON  STATION 
ALEXANDRIA,  VIRGINA 


Department  of  Defense  contractors  must  be  established 
for  DDC  services  or  have  their  ' need-to-know'  certified 
by  the  cognizant  military  agency  of  their  project  or 
contract.  All  other  persons  and  organizations  should 
apply  to  the: 


U.  S.  DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TECHNICAL  SERVICES 
WASHINGTON  25,  D.  C. 


ABSTRACT 

(Securify  class ///ca tlon  ol  title t  body  of  abstract  and  indartnjl  annotation  moat  be  entered  when  the  overall  report  la  clmaalflad) 


1  originating  activity  Appiie^  Research  Laboratory  report  security  classification 

Sylvania  Electronic  Systems,  A  Division  of  _ Unclassified _ 

Sylvania  Electric  Prod,  Inc.,  40  Sylvan  Rd,  2h-  group  <r  °t  astia  only; 

Waltham  54  .  J4as.sachuiieid:s _ _ _ _ 1 _ 

3.  REPORT  TITLE 

Admittance  of  a  Waveguide  Radiating  into  Stratified  Plasma 


4.  OESCRIPY'VE  NO.TES  (Type  ol  repost  rtno  Inclusive  dates) 

Scientific  Report  Sic ,  2 


<  name,  first  nams.  iniliit) 

inis 

ate  Auctj. s t  1 963 

7.  TOTAL  NO  OF  PAGES  59 

F.PORT  ;jo!SI 

3a.  COM  TRACT  OR  GRANT  NO 

AF1S(620; -2410 

b.  PROFIT  MO  '  £  A  "J 

c-  -atkno  464202 

(Any  olho'f  mt tft&l  me v  be 

„  94 

15.  SUPPLEMENTARY  MOTES  (For  ASTIA  use  only) 

MTS  (F or  AS 


KEY  TERMS  -  UNCLASSIFIED  ONLY 


Reentry  Commun i ce t icr 

Layers  of  Lossy 
Dielectric _ 

; 

Stratified  Plasma 


A  DESCRIPTORS  (For  ASTIA  use  only) 


15.  IDENTIFIERS  *  UNCLASSIFIED  ONLY  (t>  Model  numbers;  weapon  system ,  project ,  chemical  compound  and  trade  names 


ASTIA  FORM  28  (Test) 

FEQ  63 


6  BOOY  OF  A85TKAC”  (Including  1 1  applicable ,  f-ni pose;  M4tl>od  or  Approach;  Results;  Conclusions,  applications  and 
recommendations.) 

A  slot  covered  by  a  stratified  plasma  is  assumed  to  radiate 
into  a  wide  waveguide  instead  of  free  space.  The  slot  admittance 
approximates  the  free  space  admittance  of  the  slot  for  waveguide 
diameters  exceeding  6  to  .10  \ .  For  thick  plasma  layers  the  com¬ 
uted  slot  admittance  checks  with  earlier  admittance  calculations 
or  a  laterally  unbounded  plasma.  When  approximating  a  plasma 
profile  of  a  typical  hypersonic  re-entry,  a  multi-layer  plasma 
model  in  a  wide  waveguide  appears  to  provide  a  more  accurate  slot 
admittance  than  a  single  layer  approximation  in  a  laterally  un¬ 
bounded  geometry. 


17.  INCFXINO  ANNOTATION 

The  admittance  is  obtained  for  a  rectangular  slot  antenna 
covered  by  a  stratified  plasma  approximating  a  typical  re-entry, 
radiating  into  a  wide  rectangular  waveguide  instead  of  free  space. 


ADMITTANCE  OP  A  WAVEGUIDE  RADIATING 
INTO  STRATIFIED  PLASMA 


I .  INTRODUCTION 

The  admittance  of  a  vaveguide  or  cavity  backed  slot  which  radiates  into  a 
homogeneous  isotropic  plasma  or  dielectric  layer  of  a  thickness  that  is  larger 
than  the  wavelength  has  been  recently  computed  by  Gale.js  [l],  This  paper  also 
contains  numerous  references  to  earlier  work  on  related  problems.  The  formulation 
[l]  is  rather  involved  and  it  is  not  suited  for  determining  the  admittance  for 
thin  or  stratified  plasma  layers. 

In  this  latter  problem  it  may  be  advantageous  to  invoke  geometrical  approxi¬ 
mations  and  to  consider  the  admittance  of  a  plasma  or  dielectric  covered  slot 
which  radiates  into  a  wide  waveguide  instead  of  free  space.  Cohn  and  Flesher  [2] 
have  discussed  radiation  from  a  coaxial  line  into  a  waveguide  the  walls  of  which 
are  removed  to  infinity,  and  their  results  are  in  agreement  with  other  analyses 
not  employing  such  an  artificial  boundary  [Levine  and  Papas,  3] •  The  .waveguide 
boundary  greatly  simplifies  the  admittance  calculations  for  a  plasma  or  dielectric 
covered  slot,  if  it  is  permissible  to  use  a  waveguide  of  finite  diameter. 

Heuristic  reasoning  suggests  that  the  fields  in  the  vicinity  of  the  slot  should 
not  be  affected  by  the  presence  of  the  waveguide  walls  if  the  distance  between 
the  walls  is  large  enough.  Hence,  this  method  can  be  expected  to  yield  realistic 
susceptances  and  also  realistic  estimates  of  losses  in  the  plasma  layer.  However, 
the  radiation  fields  of  the  slot  will  be  modified  by  the  waveguide  walls.  This 
model  is  therefore  not  suitable  for  determining  radiation  patterns  of  the  slot 
and  criteria  must  be  developed  for  establishing  the  validity  of  conductance  esti¬ 
mates.  Also  a  large  number  of  terms  is  required  in  the  field  representations  for 
large  waveguide  cross-sections,  which  makes  this  approach  useful  only  if  a  high- 
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speed  computer  is  available  for  the  numerical  work. 

The  specific  problem  to  be  considered  in  this  paper  is  the  admittance  of 
a  rectangular  waveguide ,  excited  in  its  principal  mode  which  radiates  into  a 
wide  waveguide  containing  layered  isotropic  plasma  of  arbitrary  parameters  -  For 
the  purposes  of  this  paper  the  plasma  is  represented  as  a  lossy  dielectric  of  a 
relative  dielectric  constant  less  than  unity.  This  simple  representation  is 
valid  only  in  low  power  applications .  The  analytical  formulation  of  the  slot 
admittance  can  be  considered  as  straight  forward,  and  only  the  principal  results 
of  the  development  are  summarized  in  Section  II.  Convergence  considerations  and 
the  numerical  results  will  be  discussed  in  Section  III,  which  contains  also  a 
detailed  comparison  with  admittance  calculations  for  a  laterally  unbounded  plasma 
layer  [1]. 

II.  SLOT  ADMITTANCE 

2.1  The  general  formulation 

The  waveguide  geometry  under  consideration  is  depicted  in  Figure  1. 

A  small  rectangular  waveguide  of  cross-sectional  area  x_^y_^  is  joined  symmetri¬ 
cally  to  a  larger  waveguide  of  area  xoyQ>  by  a  symmetrical  aperture  of  width  2e 
and  length  26 .  The  smaller  guide  is  excited  in  the  principal  (TE  q)  mode,  but  the 
change  of  the  waveguide  cross-section  and  the  dielectric  discontinuities  generate 
a  large  number  of  TE  and  TM  modes  which  must  be  considered  in  the  formulation  of 
the  slot  admittance .  Such  calculations  have  been  carried  out  in  Appendix  A  and 
the  following  stationary  expression  for  the  normalized  slot  admittance  y^  is 


derived: 
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E  and  E  refer  to  the  aperture  fields  at  z  =  0  and  j  =  i  or  p. 
x  y  x  ^ 

coordinates  are  used  with  subscript  i,  and  the  unprimed  with  p. 
are  defined  as 
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The  prime  on  the  double  summation  (l)  designates  the  omission  of  its  n  =  1, 

m  =  0  term.  The  reflection  coefficients  of  the  TE  modes  R  and  of  the  TM 

a 

modes  depend  on  the  dielectric  structure  of  the  larger  waveguide.  For 

stratified  dielectric  layers  shown  in  Figure  2  the  reflection  coefficient  R  . 

(d  =  a  or  b)  in  the  region  of  zj+1  <  z  <  z^  is  related  to  the  reflection 

coefficient  Rn/  ,  of  the  region  z.  <  z  <  z,  .  by  the  expressions 
dlJ-1)  J  -  -  J-l 


(6) 
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R^  is  derived  by  considering  the  scalar  functions  (5 4)  and  (55)  in  two  adjacent 
layers  and  by  requiring  the  tangential  field  components  (56),  (57),  (59)  and 
(60)  to  be  continuous  across  the  interface.  The  computations  should  start  with 
j  =  1  where  =  =  0  in  (5)  and  (6).  A  series  of  computations  gives 

then  R^,  R^,  R^„  and  finally  R^  as  R&  or  in  (2)  and  (3). 

This  completes  the  formal  specification  of  the  slot  admittance  y  in 

s 

terms  of  the  electric  field  components  E^  and  E^  in  the  aperture  at  z  =0,  and 
it  remains  to  find  a  compatible  set  of  E_^  and  E^  for  the  aperture. 

In  deriving  the  slot  admittance  the  continuity  of  the  tangential  electric 
field  components  is  assured  by  the  definition  of  mode  amplitudes  following  (65) 
and  (66).  The  continuity  of  the  magnetic  field  components  and  is  specified 
by  (68)  and  (69).  The  above  equations  still  apply  to  the  derivation  of  y  as  E 

S  X 

is  made  very  small,  and  the  amplitude  of  the  TE  and  TM  modes  of  the  waveguide  may 


be  selected  in  such  a  way  that 


ff. 


E 


is  negligible  relative  to 


UfiJ 


in  (l).  The  admittance  is  now  determined  by  specifying  only  a  single  electric 
field  component  E^  in  the  aperture,  which  is  in  line  with  earlier  work  on  radia¬ 
tion  from  waveguide  backed  slots  in  absence  of  dielectric  discontinuities 
(Lewin  [4],  Stevenson  [5],  Watson  [6]).  E^  =  0  has  been  also  used  in  computations 
for  radiation  across  a  dielectric  discontinuity  from  an  oil  filled  guide  into 
free  space,  and  the  accuracy  of  such  calculations  have  been  checked  experimentally 
(Figures  4  and  5  of  Cohen,  Crowley  and  Levis  [7]).  It  may  be  noted  that  the 


continuity  condition  of  H  in  (69)  leads  after  multiplication  with  E  and  an 

y  /  f  p  \2 

(  '  'E  )  in  (l).  The  continuity  of  H  is 


integration  to  the  terms  proportional 


not  required  for  deriving  y  with  E  =0,  and  it  is  possible  that  simple  trial 

S  X 

functions  of  the  aperture  fields  which  violate  the  condition  of  H  continuity 

/  rr  \2  y 

still  lead  to  nearly  correct  results  of  y  ,  if  '  E__ ,  is  equal  to  zero  in 

®  '  1. '  i.  x  j 
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this  approximation.*  The  problem  of  non-physical  solutions  is  avoided  if  a 
neglect  of  (  f  /e  is  interpreted  as  E  negligible,  but  non-zero  in  the 
cases  where  =  0  is  not  a  strict  solution.  The  calculations  of  Section  2.2 
are  based  on  (l)  with  E  negligible  relative  to  E  . 

x  y 

It  is  also  possible  to  select  reasonable  trial  functions  with  E  £  0. 

x 

When  specifying  the  aperture  field  component  as  that  of  the  principal  wave¬ 
guide  mode  (Hx  =  Hq  sin  @xilx' )  in  a  completely  open  smaller  guide,  E^  £  0  if 
y  variations  are  allowed  for  the  E^  component.  In  this  particular  case  the 
E  integrals  of  the  slot  admittance  in  (l)  may  be  expressed  in  terms  of  E 

x  y 

integrals,  and  a  separate  trial  function  is  not  required  for  the  E^  component. 

This  development,  which  is  carried  out  in  Appendix  B,  is  summarized  in 

Section  2.3.  It  is  shown  that  the  presence  of  the  E^  component  gives  only  a 

minor  correction  to  the  real  part  of  y  ,  but  the  net  change  of  y  remains  to 

s  s 

be  determined  by  numerical  computations. 

2.2  Slot  Admittance  with  a  Negligible  E  component. 

X 

With  E  (x,y)  negligible  relative  to  E  (x,y)  the  electric  fields 
x  y 

in  the  slot  are  assumed  to  be  of  the  form 


Ey(x,y)  = 


AfA(x)  +  BfB(x) 


f(y) 


(7) 


where  A  and  B  are  the  undetermined  amplitudes  of  the  trial  functions  f.(x)  and 

A 

f  (x)  and  where  f(y)  is  yet  arbitrary.  Substitution  of  (7)  in  (l)  results 
£> 


* 

Such  an  approximate  solution  is  obtained  by  assuming  E^  =  0  in  a  plasma  layer  of 

infinite  thickness,  and  an  aperture  of  of  28  =  x_^,  2e  =  y^  when  E^  =  f(y)  violates 

the  condition  of  continuous  (See  Appendix  D) .  However,  numerical  calculations 

show  that  y  is  nearly  the  same  in  this  case  as  for  E  =  0,  E  £  f(y)  and  E  £  0. 

J  s  x  7  y  w  x  r  ’ 

E  =  f(y)  where  H  is  continuous,  (see  Figure  15). 

y  y 
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and  where  G.  and  G  represent  the  two  double  summations  of  (70). 
i  ]? 

Considering  the  ratio  of  (a/b)  as  a  variable  in  (12),  the  requirement  that  y  be 

s 


stationary  leads  to 
dy 


d( A/B)  ° 


(12) 


Equation  (12)  is  solved  for  A/B 


A  =  7BB  FA  ~  7ABFB 

B  =  7aafb-7abfa 


(13) 


Substituting  (13)  in  (8)  it  follows  that 
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with 
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The  slot  admittance  (ll)  has  been  computed  for  the  field  distribution  E  (x,y) 

y 

of  (7)  with  "optimized"  coefficients  A  and  B.  The  trial  functions  of  (7) 
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should  satisfy  the  condition 


x 


+  l) 
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Suitable  trial  functions  are  [1,8] 


and 


(16) 


(17) 


(18) 


The  parameter  k  of  the  trial  functions  in  (17)  and  (l8)  may  be  made  equal  to 
k^,  k^,  kQ  or  kept  arbitrary.  The  y  variation  of  E^(x,y)  is  neglected  and  f(y) 
becomes 

f(y)  =  -  ^  (19) 

Evaluation  of  the  integrals  (10)  and  (ll)  results  in 
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In  the  n  and.  m  summations  n  is  restricted,  to  odd.  and  m  to  even  integers.  The 
computations  are  particularly  simple  when  the  width  of  the  slot  2e  is  equal  to 
the  height  of  the  small  guide  y  ;  of  (23)  contains  only  the  m=o  term  in  this 
case . 

The  slot  admittance  can  be  also  computed  for  the  sinusoidal  field  distri¬ 
bution  fA(x)  of  (17)  with  E  =  0  in  (7).  This  gives 
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2.3  Slot  Admittance  with  a  Non-Negligible  B  Component 

The  discussion  of  this  Section  is  restricted  to  small  guides  of 

y.  =  2e  and  x.  =  2£  with  an  assumed  magnetic  field  distribution  H  =  H  sin  p  ...x' 
•’l  i  ^  x  o  xil 

in  the  aperture  at  z  =0.  The  absence  of  higher  order  modes  of  in  the  aperture 

establishes  a  relation  between  the  amplitudes  a  and  b  of  TE  and  TM  modes 

nm  nm 

following  (52),  (53)  and  (59)*  Applying  this  result  to  (65)  and  (66)  //  E  ^ 

is  related  to  //  E  .  by 
yi 
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/;  E  .  =  -I1--- . -  //  E  .  (34) 

X1  PxiPyi  yi 

for  m  f  0.  Also  the  amplitudes  of  the  TE^  modes  of  the  smaller  guide  are  equal 

to  zero  and  E  exhibits  the  same  x-variation  as  H  but  it  may  have  a  different 
y  x 

y-variation.  With  the  relative  amplitudes  of  TE  and  TM  modes  determined,  E^  of 
(56)  can  be  also  expressed  in  terms  of  (//  E^).  Hence  //  E^  of  (l)  can  be 
related  to  E^.  These  calculations  have  been  carried  out  in  the  Appendix  C.  The 
integral  / /  is  computed  as 


2x.3x  V  n(p  2  -  k.2)(-l)  2 
100  xil  1 


2  2  2 
x  -  n  x. 
o  1 


cos  p  H 
xp 


sin  P  e 

-  -  Q(P  c) 

V  ™  -I 


where  V  is  the  potential  across  the  slot  and  Q(w)  is  defined  as 
o 

2  € 

Q(w)  =  Q(Pe)  =  -  y  /  E  (u)  cos  pu  du 

o  o  y 


u  =  y'  -  e 
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The  integrals  of  E 


yp 


,  E  .  and  D  are  expressed  as 
yi 


n-1 


If  E. 


m 

\2 


yp 


2x. x  V  (-1)  2 

«•  v  °-<v> 

O  1 


(38) 


JJ  E  .  = 
yi 


V  x.  |+1 

-tj—  (-1)  Q(^)  for  n  =  1,  m  £  o 


0 


for  n  f  1 


(39) 


D  = 


V  x./2 
o  1 


(to) 


The  above  integrals  depend  on  the  function  Q(w)  or  on  the  E^  variation  in 
the  y  direction.  The  static  solution  for  the  waveguide  profile  of  Figure  1  with 
xq  »  x^  ,  carried  out  in  Appendix  B,  shows  that  in  the  vicinity  of  the  waveguide 
edge  at  u  =  e  -  &(6  -»  o)  the  electric  field  component  E  is  represented  by 

y 


E  =  -  V  4^- 
y  o  4 


2 

3«€  6 


1/3 


(41) 


On  the  symmetry  line  of  the  guide  at  u  =  o  it  is 

V  af  °-833 

This  field  variation  will  be  approximated  by 


(42) 


E  (u)  = 

y 


+  E  4l.y  (-v  U1'3) 

2 2vl/3  '  o  ’ 


(43) 


Substituting  (43)  in  (36),  Q(w)  can  be  evaluated  using  the  integrals  (1.3.8)  and 
(1.1.5)  of  Erdelyi  [9].  A  calculation  shows  that 

Q(w)  =  A  Q^w)  +  B  Q^(w)  (44) 
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where 


QjN)  = 


21/6  r( 2/3)^ 


QgCw)  =  21//6  r(2/3)  >/^T 


w 

r  Jl/6(») 

i7T 

W  ' 


k  J7/6(v) 

7/6 


w 


(45) 


(46) 


and  Jn(w)  is  the  Bessel  function  of  first  kind  of  order  n  and  P(z)  is  the  gamma 
function.  The  constants  A  and  B  are  determined  by  requiring  that  E^  of  (4-3)  is 
equal  to  (4l)  near  the  edge  and  that  the  integral  of  E^(u)  across  the  aperture  is 
equal  to-V  (or  Q(w)  =  1  for  p  =  o) .  It  follows  that 


A  =  1  r(  7/6) 

h  r(2/3) 


21  fjJ 3  \  1 ^ 

5^  2n 


(47) 


B  =  -  A  +  | 


y?  \V3 

2rt 


(48) 


With  A  and  B  determined  as  above,  the  E  (o)  of  (43)  is 

y 


E  (o)  = 

y 


0.852 


(49) 


which  is  close  to  the  correct  static  value  of  (42).  The  slot  admittance  y  may  be 

s 

computed  after  substituting  (34)  to  (40)  into  (l).  Die  summations  are  extended 
over  even  integers  m  and  odd  integers  n,  except  forthe  primed  summation  where 
n  =  1  because  of  (39) • 

The  static  solution  for  the  slot  in  an  infinitely  thin  screen  gives 
2  2  _1/2 

Ey(u)  =  -  (Vq/j t)(e  -  u  )  and  Q(w)  =  J_  ;(w) .  This  will  give  higher  values  of 
integrals  than  the  more  accurate  forms  (43)  and  (44).  When  ignoring  the  y  varia¬ 
tion  of  E  (x,y)  in  the  aperture  Q(w)  becomes  sin  w/w  and  the  integrals  //  E_  , 
y  xp 

//  E  .and  //  E  .  are  equal  to  zero  according  to  (35),  (39)  and  (34).  The  slot 
yx  xi 
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admittance  computed  from  (l)  using  this  particular  form  of  Q(w)  in  (38)  is 
the  same  as  (33)  in  the  limit  of  k  -*■  0^^  • 

The  error  of  y  which  results  from  the  neglect  of  E  in  the  aperture 

S  X 

(E^  may  he  non-negligible  for  E^(u)  of  (43))  can  be  estimated  by  computing  the 

ratio  between  terms  proportional  to  f  J  E  .  and  //  E  . 

xj  y  j 

Using  the  small,  argument  approximations  of  sin  w/w  and  of  Q(w)  of  (4-4), 
where  w  =  0  e  ,  it  follows  that 

yp 


//  E  n[ it2  -  (k.x.  )2] 

J  J  xp  _  _ _ v  1  1' 


II  E. 


12  m 


0.032w 


(50) 


yp 


*=  0.008«2  nm  (  — 


y. 


for  w  <  1 


For  yQ/yi  >  30  ,  ff  E ff  E^  <  0.1  for  nm  <  1140.  Hiis  range  of  nm  encompasses 

the  more  significant  terms  of  (l)  and  the  terms  proportional  to  (  //  E  )  will 

xp 

cause  only  a  minor  change  in  y  ,  as  long  as  F(R)  and  G(R)  remain  of  approximately 

s 

the  same  magnitude. 

2  2 

The  ratio  of  the  terms  proportional  to  ( ff E^ )  and  (// E  )  in  the  primed 
summation  of  (l)  is  computed  with  the  aid  of  (34)  as 


k.2  -  0  -2  (Jf E  .f 

_ 1 _  yi  J  •  xi 

0  2  -  k.2  (ffE  .f 

xil  1  yx 


2  2 
k.  -  0 
1  xil 


P 


xil 


(51) 


For  x.  >  0.6l  the  above  ratio  is  more  than  0.44,  and  the  terms 
1  - 

(// E  . )  are  significant,  particularly  for  larger  values  of  x. . 

XI  1 

is  real  and  therefore  contributes  only  to  the  imaginary  part  of 
imaginary) . 


proportional  to 
This  summation 
ys  in  (1)  ( 7±1 


is 
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The  presence  of  an  E^  component  will  give  only  a  minor  correction  to 

Re  y  i  the  change  in  Im  y  will  depend  on  the  ratio  of  the  primed  summations 
s  s 

to  the  other  summations  and  will  be  examined  by  numerical  calculations. 


Ill  DISCUSSION 


The  slot  admittance  y  has  been  derived  in  Section  2.2  with  a  two  term 

5 

trial  function  for  the  aperture  field  component  E  ,  which  was  assumed  to  exhibit 
no  y -variations .  E  was  neglected  in  these  calculations.  The  numerical  results 
shown  in  Figures  3  to  11  and  13  to  14-  are  obtained  using  this  formulation.  The 
derivation  of  Section  2.3  assumes  the  principal  mode  H  distribution  in  the  slot 
and  it  allows  y-variations  of  the  E^  component.  E^  f  0  in  this  approach.  The 
results  obtained  by  using  the  two  different  formulations  are  compared  in  Table  1 
and  in  Figures  15  to  17* 

The  numerical  calculations  will  illustrate  first  the  variations  of  the 

slot  admittance  when  a  slot  radiates  into  an  open  waveguide  of  increasing  size. 

Calculations  were  made  for  xq/x^  increasing  in  increments  of  0.1  and  y^  =  xq. 

The  slot  admittance  of  Figure  3  exhibits  systematic  variations  for  smaller  values 

of  x  /x.,  but  the  admittance  points  appear  to  be  scattered  with  smaller  deviations 
o  1 

for  xq/x.  >  20,  although  there  are  several  widely  scattered  points.  Also,  there 
are  deviations  in  both  admittance  components  G  and  B,  which  is  contrary  to  the 
elementary  reasoning  presented  in  Section  I.  Examination  of  (l4)  shows  that  y 


may 


become  infinite  if  H  =  00  in  (24)  to  (27)  which  occurs  in  absence  of  plasma 


(R  =  R^  =  o)  at  the  mode  cut-off  when  7^  =  7^  -*  0.  (With  a  plasma  layer  R&,  R^  f 

If  7  i  0  and  7  -*  0  ,  R  and  R,  -*  -1,  but  the  ratios  7  /(l  +  R  )  and  (l  +  R,  )/y 
o  p  'a  b  7  p  a  d  p 

and  also  H  remain  finite.  With  a  single  plasms,  layer  and  7  -*  0,  H  may  become 
n  o  n 

infinite  only  if  exp(27  z  )  =  +  l) .  Considering  the  free  space  admittance,  H 


P  o 


rr  03 


\  l  2  2 

occurs  for  =  —  s  n  +  m  where  n  =  odd  and  m  =  even  integers .  These  values 

of  are  indicated  by  vertical  lines  in  Figure  3  for  the  smaller  values  of 

xo/x^.  The  above  values  of  x^  correspond  to  the  apparent  discontinuity  points 

of  the  slot  admittance.  The  widely  scattered  admittance  points  for  larger 

values  of  x  /x  are  also  very  close  to  these  cut-off  values  of  x  .  Although 
o  i  o 

the  separation  between  the  cut-off  points  of  the  modes  is  decreased  for  increas- 

J 

ing  values  of  x^  in  Figure  3,  one  has  to  approach  the  cut-off  point  more  closely 
to  obtain  a  sizeable  deviation  from  the  mean  of  the  calculated  admittance  figures . 
This  accounts  for  the  smaller  number  of  widely  scattered  admittance  points  for 
xq/x.  large.  The  free  space  admittance  of  the  present  slot  is  y  =  0.63  -  i  0.35 
following  Lewin  [4]  or  y  =0.65  -  i  0.38  following  Galejs  [l].  The  waveguide 
model  is  therefore  effective  for  obtaining  an  approximation  to  the  free  space 
slot  admittance  provided  the  ratio  of  slot  length  to  the  outer  waveguide 
dimension  exceeds  10.  In  order  to  exclude  the  widely  scattered  admittance  points 
from  consideration,  the  slot  admittance  is  required  to  be  a  continuous,  gradually 
changing  function  of  the  guide  dimensions  in  the  vicinity  of  the  selected  outer 
guide  dimension. 

Examples  of  admittance  calculations  are  shown  for  XD/X^  =  10.5  and  19.5 

in  Figures  4  to  9  for  a  single  layer  plasma  of  relative  dielectric  constant 

e  /e  .less  tangent  tan6  and  thickness  h  =  z  .  The  slot  admittance  starts  out 
P  o  p  P 

with  a  free  space  admittance  for  h/x  much  less  than  1  and  approximates  the  slot 
admittance  of  a  semi -infinite  plasma  layer  for  h/x  much  larger  than  1.  Although 
the  admittance  variation  is  relatively  smooth  for  h/x  much  less  than  1,  it 
becomes  oscillatory  for  h/x  larger  than  1.  For  relatively  thin  layers  of  plasma 
the  admittance  is  less  dependent  of  guide  dimensions,  if  the  plasma  layer  provides 
a  sharp  discontinuity  with  respect  to  free  space  (e.g.,  £p/£0  «  f) •  In  this 

case  the  guide  admittance  changes  are  determined  principally  by  the  effects  of 
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the  nearby  plasma  and  free  space  interface,  with  guide  walls  being  of  only  a 

minor  significance.  However,  for  thick  plasma  layers  the  plasma  and  air 

interface  is  at  a  distance  comparable  to  the  distance  from  slot  to  the  guide 

walls  and  the  admittance  depends  more  on  guide  dimensions.  Also,  there  are 

oscillations  of  the  admittance  for  h/x  larger  than  1.  The  slot  admittance  for 

a  semi-infinite  homogeneous  plasma  of  Galejs  [l],  which  is  indicated  dashed  in 

Figures  4  to  9>  is  approximately  equal  to  the  average  of  the  oscillatory 

conductance  curve  for  h  and  XQ/Xi  large.  The  admittance  curves  of  Figures  8  and  9 

for  x  =  y  =  19.5  xj  appear  to  approximate  the  admittance  of  the  semi-infinite 
00  i 

plasma  closer  than  Figures  6  and  7  for  x^  =  yQ  =10.5  x^.  For  h  >  X  the  slot 
admittance  is  about  the  same  as  the  slot  admittance  for  a  semi -infinite  plasma 
layer,  provided  the  conductivity  oscillation  is  smoothed  out.  The  slot  admittance 
tends  toward  the  free  slot  admittance  for  h  <  X/2.  The  above  conclusions  support 
the  findings  of  plasma  layers  analysis  of  Galejs  [l]  which  was  strictly  applicable 
only  to  layers  of  h  »  X. 

The  variational  approximation  to  the  voltage  distribution  along  the  slots 

is  computed  from  (7)  for  the  slot  and  plasma  parameters  of  Figures  4,  5>  8,  9  and 

are  shown  in  Figure  10.  The  voltage  distribution  is  seen  to  change  only  slightly 

with  varying  plasma  parameters.  These  calculations  are  made  with  k  =  kQ  in  the 

trial  functions  (17)  and  (l8) .  For  k  =  k  the  voltage  distribution  is  changed 

P 

for  the  same  plasma  parameters  to  the  form  indicated  in  Figure  11.  These  distri¬ 
butions  are  about  the  same  as  for  the  laterally  unbounded  plasma  in  Figure  7  of 
Galejs  [1],  The  calculated  admittance  is  nearly  the  same  for  both  values  of  k  of 
the  trial  functions,  as  may  be  seen  from  the  examples  of  Table  1. 

The  accuracy  of  the  voltage  distribution  shown  in  Figures  10  and  11  can 
be  strictly  determined  only  by  comparing  these  variational  approximations  with 
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a  more  accurate  result,  if  this  were  available.  However,  the  assumed  form  of 

trial  function  is  general  enough  to  allow  for  wide  variations  of  the  field 

distribution  along  the  slot,  such  as  observed  in  cavity  backed  slot  antennas 

[Galejs  1,  7] .  The  close  resemblance  between  the  voltage  distributions  computed 

with  k  =  k^  and  k^  suggests  that  the  true  distribution  may  be  intermediate. 

These  distributions  may  be  approximated  by  the  principal  waveguide  mode 

(E  ~  sin  jtx'  /x.  )  for  x.  =  2i  . 
y  1  1 

It  was  also  pointed  out  that  the  present  model  is  suited  for  computing 
the  effects  of  a  stratified  dielectric  or  plasma  layer.  This  will  be  shown  for 
a  re-entry  plasma  sheath  of  Figure  20  by  Rotman  and  Meltz  [8]  which  is  reproduced 
in  Figure  12.  The  relative  dielectric  constant  and  the  loss  tangent  are  computed 
from  standard  formulas  for  various  frequencies  below  and  above  the  gyrofrequency 
of  plasma.  The  slot  dimensions  are  kept  constant  in  relation  to  the  wavelength. 
The  perturbation  of  the  plasma  profile  by  the  antenna  fields  is  ignored  and  the 
slot  admittance  is  computed  for  a  seven  (7)  layer  approximation  to  the  profile  of 
Figure  12,  which  is  depicted  in  Figure  13.  The  voltage  distribution  is  shown  in 
Figure  14.  The  calculations  made  for  estimated  average  plasma  parameters 

Q  O  g 

(Ng  =2.5  x  l(x  cm  and  v  =  3-8  x  10  sec  )  are  shown  dashed  in  Figure  13. 

These  average  data  agree  with  computations  made  for  a  laterally  unbounded  plasma 
layer  [l]  but  they  provide  significant  errors  relative  to  the  multilayer  model 
particularly  for  frequencies  below  the  gyrofrequency  of  the  plasma. 

The  preceeding  numerical  calculations  were  based  on  electric  aperture 


fields  with  neglected  E  and  with  a  y-independent  E  .  These  calculations  are 

x  y 

compared  in  Table  1  with  calculations  made  for  the  principal  mode  H  field  in  the 

aperture,  where  E  exhibits  the  static  y-variation  and  E  ^  0.  The  differences 
’  y  x 

between  the  two  sets  of  calculations  are  small.  Even  in  the  case  where  E  k  0, 

x  ’ 
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the  neglect  of  it  does  not  affect  the  slot  conductance  and  changes  the  susceptance 
by  less  than  1$.  The  change  from  the  static  E  (y)  to  a  y-independent  E  causes 

y  y 

a  negligible  change  in  G  and  about  7$  change  in  B. 

Tlie  changes  of  the  susceptance  B  are  larger  than  the  changes  of  G  for  the 

cases  listed  in  Table  1  and  further  susceptance  calculations  are  shown  in 

Figure  15  for  the  same  plasma  and  antenna  parameters  as  in  Figure  7-  There  is  a 

change  of  B  caused  by  a  change  of  E  (y),  but  none  by  the  neglect  of  E  within 

y  x 

the  accuracy  of  the  plot. 

The  slot  admittance  of  an  open  waveguide  (x^  =  2i>  ,  y^  =  2e)  radiating  into 

a  laterally  unbounded  plasma  has  been  recently  calculated  by  Villeneuve  [ll].  He 

assumes  the  principal  waveguide  mode  as  the  aperture  field  E^  with  E^  =  0,  and  the 

slot  admittance  which  is  obtained  after  numerically  evaluating  a  two-fold  Fourier 

transform  should  correspond  to  the  admittance  computed  from  (33)  in  the  limit  of 

k  -*  Pxil  =  rt/Xjy  A  comparison  of  the  results  by  Villeneuve  [ll]  with  computations 

following  Section  2.2  (E  f  f (y),  E  neglected)  and  of  Section  2.3  (H  of  the 

y  x  x 

principal  waveguide  mode  in  the  aperture,  static  E^(y)  and  E^  f  0  or  =  0)  is  shown 

in  Figures  1 6  and  17.  The  free  space  slot  admittance  of  Lewis  [4]  or  Cohen  et  al. 

[7]  is  also  indicated  in  Figures  1 6  and  17.  The  results  of  Villeneuve  [ll]  agree 

with  those  for  E  £  f(y)  of  the  smaller  values  of  e  ,  but  there  are  differences 

y  P 

up  to  5 1°  for  e  /e  >0.6.  The  neglect  of  E  in  the  solution  with  the  static 
r  '  P  o  x 

E  (y)  dependence  gives  a  detectable  change  in  the  susceptance  of  Figure  17  because 
of  the  increase  in  x^/l  relative  to  Figure  15.  (This  is  in  agreement  with  the 
discussion  following  (51)). 

The  comparison  of  slot  admittances  in  Table  1  or  Figures  15  to  17  which 
are  based  on  principal  mode  electric  or  magnetic  fields  as  the  aperture  trial 
functions  may  be  also  of  interest  in  connection  with  a  recent  discussion  of  these 
Transactions  (Unz  and  Hodara,  [9]). 
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The  above  work  has  shown  that  the  plasma  layers  in  a  large  waveguide 

may  be  considered  as  a  substitute  for  laterally  unbounded  plasma  layers.  The 

waveguide  model  makes  the  problem  mathematically  simple  and  suitable  for 

computer  calculation.  Hie  dimensions  of  the  large  guide  may  be  selected  to  be 

such  that  the  slot  has  approximately  the  correct  admittance  in  the  absence  of 

plasma  layers.  The  slot  will  exhibit  the  correct  admittance  behavior  of  a  slot 

< 

in  an  infinite  screen  for  thin  plasma  layers  (h  <  X) .  Furthermore,  the  dimensions 
of  the  large  guide  are  less  critical  if  the  electrical  characteristics  of  the  thin 
plasma  layer  differ  significantly  from  those  of  free  space  (e^/e^  «  1  or 
tan  8  »  l).  For  thick  plasma  layers  (h  >  X)  the  calculated  conductance  exhibits 
oscillations  which  are  absent  in  laterally  unbounded  plasma  layers.  However,  the 
average  of  this  oscillatory  plasma  conductance  gives  the  correct  admittance 
figure  in  the  cases  discussed  earlier.  It  appears  that  a  multilayer  plasma  model 
in  a  wide  waveguide  provides  a  more  accurate  slot  admittance  than  a  single  layer 
approximation  of  a  plasma  profile  in  a  laterally  unbounded  geometry.  Most  of  the 
calculations  make  use  of  a  two  term  trial  function  for  the  slot  fields,  but  the 
fields  do  not  differ  significantly  from  the  principal  waveguide  mode  and  the  trial 
function  and  hence  the  calculations  may  be  further  simplified. 

The  technique  described  in  this  paper  is  applicable  for  obtaining  admit¬ 
tance  estimates  for  slot  antennas  in  presence  of  layered  isotropic  plasmas  or 
dielectrics.  At  present,  there  are  no  such  treatments  for  a  laterally  unbounded 
plasma.  The  present  calculations  will  be  correlated  with  antenna  impedance 
measurements  on  hypersonic  reentry  vehicles  and  also  with  laboratory  measurements 
on  a  glow  discharge.  The  multi-layer  formalism  accounts  also  for  glass  wall 
effects  of  discharge  vessels.  The  lateral  bounds  of  the  plasma  layer  assumed  in 
the  analysis  appear  of  no  consequence  in  such  applications.  The  waveguide  model 
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has  been  also  applied  to  plasma  or  dielectric  covered  cavity  backed  annular 
slots  [Galejs,  13].  This  model  is  useful  for  representing  debris  covered 
hardened  slot  antennas  operating  in  the  LF  or  HF  bands. 
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APPENDIX  A 


SLOT  ADMITTANCE 

The  waveguide  geometry  under  consideration  is  depicted  in  Figure  1. 

The  smaller  guide  (z  <  0)  is  excited  by  waves  in  its  principal  TEL^  mode  which 

propagates  in  the  positive  z -direction.  There  are  also  reflected  waves  in  this 

guide.  Waves  propagating  in  the  positive  and  negative  z  directions  are  excited 

in  the  larger  guide  (z  >  0).  The  fields  in  the  guides  consist  of  a  superposition 

of  TE  and  TM  modes  that  can  be  derived  from  the  scalar  functions  I.  and  $ . 

J  J 

respectively.  The  subscripts  j=i  or  p  designate  the  smaller  guide  or  the  initial 
portion  of  the  larger  guide.  For  an  assumed  exp(-i^t)  time  dependence  of  the 
fields  it  follows  that 

7  n  z  -7 . ,  Z 

T.  =  a,„  cos  p  x'(e  +Re  ) 

1  10  xil  ' 


/  / 

n  m 


a  cos  p  .x'  cos  B  .y'  e 
nm  xi  Kya/ 


-7.  z 
i 


(52) 


<i>.  = 
1 


I 


v  r  _7iz 

b  sin  3  .x'  sin  B  .y'  e 
>  nm  xi  yaf 


n  m 


r->  7  z  -7  z 

A  cos  p  x  cos  Py(e'P+Re-P  ) 

nw  '  vr\  1  '  «  ' 


nm  xp 


yp 


(53) 


(54) 


n  m 


7  z  -7  z 

B  sin  3  x  sin  p  y(e  ^  +  Re  P  ) 

nm  xp  yp  b  ' 


n  m 


(55) 
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where 


13  .  =  22 

XI  X. 

X 


p  .,  =  — 

Xll  X . 


xp  X 


7 .  =  i  .  /  k2  -  p2  .  -  p2. 
J  M  J  xj  yj 


k .  =  .  [~  m2u  e  .  +  i^i  a . 

J  V  °  J  °  J 


rii=  \K-4i 


R,  R&  and  R^  are  reflection  coefficients.  The  prime  on  the  double  summation  (52) 
designates  the  omission  of  its  n=l,  m=0  term.  The  field  components  are  related 


to  ? .  and  $ .  by 
J  J 


E  .  =  -2-  y.  +  —  ®. 

xj  Cy  j  oz  ox  j 


e  .-  -  4-  y .  +  A-  «. 

yj  ox  j  oz  ay  j 


E  .  = 

zj 


H  .  = 
xj 


2  d 

y.  +  k .  4- 

J  J  oy  J 
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(60) 


H 


yj 


6z  <3 y 


- 

J 


k2 

J 


o 

cix 


(61) 


After  multiplying  (57)  with  sin  p  and  integrating  over  the  aperture  at 

z=0  it  follows  that 


anr  =  - - rv.'b'v  II  E  sm  3  . nx'  dx'  dy' 

10  B  ...x.y.tl+R)  JJ  y  xil 

Xll  11 


(62) 


Equation  (56)  is  multiplied  with  cos  3  .x  sin  3  .y  and  (57)  with  sin  3  -x  cos  3  ,y. 

yj  yj 

After  integrating  over  the  aperture  this  gives  two  pairs  of  equations  (j=i  or  p) 


with  the  unknowns  a  and  b  for  j=i  or  A  and  B  for  j=p.  After  introducing 

nm  nm  nm  nm 

the  abbreviations 


II  =  II  \  cos  3xjx  sin  3yjy  dx  dy  (63) 


and 

[[  E  .  =  II  E  sin  3  -x  cos  3  .y  dx  dy  (6k) 

yj  y  xj  yj 


where  E^  and  E^  refers  to  the  aperture  fields  z=C,  these  coefficients  may  be 
computed  as 


a 

nm 


b 

nm 


2e 


m 


x.y.(32.  +  32. ) 

1  1  xi  ^yi 


P  .  //  E  .  -  3  .  JJ  E  . 
yi  J  xi  xi  J J  yi 


2  e 


m 


i*  ■  II  K<  +  f*  •  II  e 


x . y . ( p2 .  +  3d. )7.  L  X1  X1  yi  yi  J 
i^iVrxi  yi'  1 


(65) 

(66) 
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where 


e 

m 


1  for  m=0 

2  for  m^O 


The  subscripts  i  should  be  replaced  by  p  for  obtaining  A  and  B  from  (65) 

nm  nm 

and  (66)  respectively.  Also,  the  right  hand  sides  of  (65)  and  (66)  should  be 

multiplied  with  (l  +  R  )  ^  and  (R^  -  l)  1  respectively.  After  substituting 

a  ,  b  ,  A  and  B  in  (52)  to  (55) »  the  magnetic  field  components  H  .  and 
nm’  rmr  nm  nm  xj 

H  .  of  (59)  and  (60)  are  related  to  the  integrals  of  the  electric  field  components 

y  j 

(63)  and  (64).  After  equating  the  tangential  magnetic  field  components  in  the 
aperture  one  obtains  two  coupled  equations.  Letting 


D  =  //  E  sin  3  ..x'  dx'  dy' 

J  J  y  xil 


(67) 


it  follows  that 


itf  d  sin  3xii?  =  H  Ey(x,y^Gi  ^  dy  +  M  Ex(x»y)G2  dy 

(68) 


II  Ey(x,y)G2  dx  dy  =  //  Ex(x,y)G3  dx  dy 


(69) 


where 


r-t  I  ’  €  ?  2 

G,  =  -  (3  .  -  k.)  sin  3  .x'  cos  3  .y'  sin  3  .6'  cos  3  .7' 

t  --  --  —  ' '  -vi  -1  -v -1  r  n  w  1  -v i  1  1 


n  m 


Xiyi7i  X1  1 


xi 


yi 


XI" 


yi 


+ 


m 


/ 

n  in 


xy(p2  +  32 
oc'^xp  yp 


7  3 


2  1  -Ra 


p  xp  1  +  R_ 


0  kc  1  +  R 

c>2  _p  _ b 

'  r.  1  -  R,,  J 


sin  S  x  cos  3  y  sin  3  t  cos  3  R 
xp  7PJ  xpb  yp 


(70) 
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n  n  ’e  3  .3  . 

)  )  — — —  sin  3  , x 1  cos  3  .y'  cos  3  -S'  sin  3  .R1 

'  Z_  xiy-?7i  xi  yl  X1  yi 


n  m 


e  3  3 
m  xp'yp 


xv. 


(3  +  3^  ) 


1  -  R 


a 


P  1  +  R, 


n  m  o'  o'  xp  yp 


a 


k;  1  +  Rb 

r  1  -  ^  j 


sin  3  x  cos  3  y  cos  3  I  sin  3  R 
xp  ypJ  xps  yp 1 


(71) 


n  m 


6  (k?  -  32.) 

m  i  yi 

x.y.7 . 
i  i  i 


cos  3  .x1  sin  3  .y’  cos  3  .1'  sin  3  . -n ’ 
xi  lyiJ  yi 1 


m 


x  y  (32  +  3^  ) 

n  m  oJo  xp  yp 


-7_  3! 


2  ^a 


P  ^yp  1  +  R 


a 


2  2 

3_  1  +  Rb 


*£  g 


1  - 


*b  J 


cos  3  x  sin  3  y  cos  3  £  sin  3  R 

xp  yp  xpb  yp 


(72) 


The  normalized  slot  admittance  of  the  waveguide  is 


1  -  R 
1  +  R 


(73) 


A  variational  principle  for  computing  y  may  be  constructed  by  multiplying 

s 

(68)  with  E  and  (69)  with  E  and  by  integrating  over  the  slot.  The  two  terms  of 
y  x 

(68)  and  (69)  proportional  to  G^  are  equal  and 


ys  = 


fill  Ey(x,y)  Ey(|,T))  G1  dx  dy  d|  dq  +  ffff  Ejx,y)  Ex(|,tj)  G3  dx  dy  dg  dq 

(74) 


7il  2 
— —  D 
x.y. 

1  1 
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The  stationarity  of  (7^)  can  be  demonstrated  by  examining  the  first 

order  change  Sy  due  to  small  changes  SE  and  SE  about  their  correct  values 
s  x  y 

E  and  E  .  As  a  result  of  these  calculations  by  is  shown  to  be  zero, 
x  y  s 

Because  of  the  symmetry  of  G^  and  the  variation  of  yg  due  to  small 

changes  SE  and  SE  can  be  written  as 
6  y  x 


fVi 

hi®2 


JJJJ  E  SE  G  +  ////  E  SE  G 
y  y  1  J  J  J  J  x  x 


JJJJ  EEG  +  JJJJ  E  E  G0 
yyl  xx3 


II  sin  PxilS  dT) 


/D 


(75) 

where  the  differentials  dx  dy  d|  dr]  have  been  omitted  from  the  four  fold  integrals 
for  the  sale  of  brevity.  Equation  (68)  is  first  multiplied  with  SE  (£,q)  and 

y 

integrated  over  the  aperture.  Then  (68)  is  multiplied  with  E  (g,Tj)  JJ  SE  (£,q) 

y  y 

sin  <3-1  and  again  integrated  over  the  aperture.  Combining  the  resulting 

two  equations  it  follows  that 


JJJJ  E  SE  G-.  +  JJJJ  E  SE  G  = 
y  yl  x  y  2 


JJJJ  EEG  +  JJJJ  EEG 
J  J  yyl  x  y  2 


II  6EyU,q)  sin  Pxil5  d|  dq 
D 


(76) 


Taking  the  variation  of  (69),  multiplying  the  resulting  equation  with  Ex(£,q) 
and  integrating  over  the  aperture  results  in 

JJJJ  ExSEyG2  =  JJJJ  ExSExG3  (77) 

Multiplying  (69)  with  Ex  (|,q)  and  integrating  over  the  aperture 

JIJJ  EEyG2  =  (73) 
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After  using  (77)  and  (78)  to  eliminate  G^  from  (7 6),  the  result  shows  that 

5y=  of  (75)  is  equal  to  zero, 
s 
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APPENDIX  B 

STATIC  SOLUTION  FOR  THE  WAVEGUIDE  PROFILE 


The  waveguide  profile  of  Figure  1  can  be  represented  by  the  t-plane 
(t  =  z  +  iy)  configuration  of  Figure  18  if  a  conducting  plane  is  substituted 
for  the  symmetry  plane  y  =  y  / 2  of  Figure  1  and  if  the  effects  of  the  outer 
guide  wall  are  neglected  (x  »  x^  ) .  The  space  between  the  electrodes  of  the 
t-plane  is  mapped  into  the  upper-half  of  the  w-plane  (w  =  u  +  iv)  by  the 
transformation 


dt 

dw 


w 


w  -  i 


(79) 


where  C,  is  a  constant.  Integrating  (79)  between  the  two  B  points  of  the  t-plane 
and  around  the  B  point  of  the  w-plane  gives  the  constant  as  =  ie/n.  Inte¬ 
grating  (79)  and  noting  that  the  t  =  ie  point  is  mapped  into  w  =  1  results  in 


The  complex  potential  function  for  the  electrodes  of  the  w-plane  is 

V 

W  [w(t)]  =  -  §«  •£n  w  (81) 

The  potential  is  the  imaginary  part  of  W  and  the  electric  field  is  given  by 
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where  the  asterisk  indicates  a  complex  conjugate.  The  electric  field  component 
Ey  can  be  determined  as  the  imaginary  part  of  (82)  in  the  z  =  0  plane  for 
0  <  iy  <  ie  only  after  inverting  the  mapping  function  (80).  The  general 
inversion  problem  will  not  be  attempted  here,  and  the  consideration  will  be 
restricted  to  points  in  the  vicinity  of  the  corner  t  =  i€  and  on  the  ground 
plane  t  =  0. 

The  point  w  =  1  corresponds  to  t  =  ie  and  s/w-1  will  be  small  for 
t  =  i(e  -  6)  with  &  -*•  0.  Expanding  the  inverse  tangent  of  (80)  results  in 


Substituting  (84)  in  (82)  gives 


E  +  i  E 

z  y 


-\)lh  ( 

3«  6  ' 


(85) 


(The  ratio  of  E^  and  iE^  depends  on  the  direction  along  which  the  corner  is 
approached.  When  approaching  along  the  symmetry  line  of  the  corner,  t  =  ie  +  &(l-i) 
and  a  similar  calculation  shows  that  Ez  +  i  E  *»  (l  -  i),  which  follows  also  from 
elementary  considerations). 


The  t  =  0  point  is  mapped  on  the  negative  real  axis  of  the  w-plane  and  the 
corresponding  value  of  w  can  be  determined  from  (80)  as  a  solution  of  the 
transcendental  equation 


-  tan"1  (i£)  =  -  | 

(86) 

or 

where 

i 

=  coth  £ 

(87) 

i£  =  ‘v/w  -  1 


-29- 


The  solution  of  (87)  can  be  found  from  tables  as  £  =  I.1997  «  1.2  or 
v  =  -  0.44.  Substituting  this  value  of  w  in  (82)  the  electric  field  is 
computed  as 


E  +  i  E 
x  y 


V 

o 

2e 


1 

i  1.2 


V 

i  ~  0,833 


(88) 
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APPENDIX  c 


FIELD  COMPONENTS  ASSOCIATED  WITH  THE  PRINCIPAL  WAVEGUIDE  MODE  AS  THE  H 

COMPONENT  IN  THE  APERTURE 


x 


The  magnetic  field  component  H  .  can  be  computed  from  (52),  (53)  and 

XI 


(59)*  Tie  aperture  field  component  H  ^  is  given  by 


tr  it  • 

H  .  =  H  sin  - 

xi  o  x. 

1 


(89) 


over  the  aperture  of  the  small  guide  (y^  =  2e,  x_^  =  22)  only  if 


a  6.7.  +  b  k.  8  .  =  0 

nm  xi  1  nm  1  1  yi 


(90) 


for  nm  /=  10.  The  coefficients  a  and  b  are  defined  by  (65)  and  (66)  in  terms 


nm 


nm 


of  E  and  E  integrals  (63)  and  (6L).  Substituting  (65)  and  (66)  into  (90)  the 
x  y 

integrals  E  ^  are  related  to  integrals  of  E^  as 


8  2  -  k.2 

Exi  -  -fry-  "  Eyi 
xi  yi  J 


(91) 


where  m  k-  0.  For  m  =  0.  £3  .  =0  and  a  =0  for  n  p  1  according  to  (90). 
'  ’  yi  no  • 

ano  of  (65)  is  zero  only  if  ff  E  =  0,  -which  implies  that 


But 


E  ,(x,y)  =  E(y)  sin  — 

y 1 


(92) 


After  substituting  (9l)  in  (65)  and  (66)  the  coefficients  a  and  b  can  be 


nm 


nm 


expressed  in  terms  of  ff  E  ^  only  and  E^  may  be  computed  from  (52),  (53)  and 
(56)  as 


E  .  = 
xi 


2  2  "N 

3  -  kf  1 

xii  1  7  r 


q=2 


x.y.  8  ...  8  • 
iJi  xil  yi 


ff  E  .  sin  3  .-.x'  cos  3  .y'  dx'  dy' 
J  J  yi  xi  1  yi 


.  cos  3  -tx'  sin  3  .y' 
xil  yi’7 


(93) 


-31- 


where  the  definition  of  p  ^  is  changed  to  qjt/y\  in  order  to  distinguish  the 

index  of  the  summation  from  the  constant  value  of  m  in  6  =  mjt/y  .  Because 

yp  o 

of  the  continuity  of  E  in  the  aperture,  the  integrals  ff  E  of  the  y 

x  xp  s 

expression  may  be  evaluated  as 


ff  E  =  //  E  .  cos  3  x  sin  8  y  dx  dy  (94) 

xp  J  J  xi  xp  ypJ  J  ' 

Equation  (93)  is  substituted  in  (9*0  and  the  x',  x  and  y  integrals  are  seen 
to  be  elementary.  The  q-summation  is  evaluated  prior  to  the  y'  integral  with 
the  aid  of  the  relation 


V  cos  2n 6 


n=l 


(2n)2-a2 


1 

2a 


1  jt  cos  a 6  +  cos  a(n -6) 

a  2  sin  an 


(95) 


which  has  been  derived  from  (558)  of  Jolley  [l4]  and  Ex.  3  p.  371  of  Bromwich 
(15].  After  carrying  out  the  above  steps  (94)  becomes 


ff  E 


8  x.\ 


n-1 


S  +  1 

2 


1  o  /_  2  .  2v  (-1)  2  n  ro  Xin  /  , 

xp  2  ^xil  ki  ^  2  2  2  cos  [|3  2^-1^  S 

y. Jt  x  -  n  x.  ^ 

1  °  1  (96) 


where 


y. 


sin(p  ypyi/2)  ,  1 


23 


yp 

3 

yp 


I  ay-  Vy,){gL--  T 

o  ^  ^  yp 


l+cos3  y. 


cos  3  y'  — r 


yp  x 


yp  singly. 


+  sin  3  y' 

yp 


(97) 


The  numerical  value  of  S  and  of  ff  E  depends  on  the  estimate  of  E  (y'). 

xp  yw 

For  E  (y' )  =  const,  the  y'  integral  of  (97)  is  equal  to  zero,  which  can  be  seen 

y 

more  readily  from  (93)*  This  makes  //  E  =0  and  the  solution  is  reduced  to  the 

xp 

form  discussed  in  Section  2.2.  The  static  solution  of  Appendix  B  for  the  wave¬ 
guide  profile  along  the  y'  =  y^/2  plane  provides  another  estimate  of  E  (y'). 
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If  E  (y* )  is  an  even  function  about  the  slot  center  (y1  =  yj2) , 
S  can  be  expressed  with  the  aid  of  Q(w)  of  (36)  as 


S  = 


Vi 

43 


yp 


sin(Q  e) 

7a - —  +  Q(3  e) 

(6ype)  Vp 


which  completes  the  evaluation  of  //  E  . 


//  E  ,  fj  E  .  and  D  are  evaluated  after  substituting  (92).  The 

yp  yi 

results  of  these  integrations  are  shown  as  (38),  (39)  and  (40). 


(98) 
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APPENDIX  D 


FIELD  COMPONENTS  FOR  E  =0 

x 


Tiie  field  components  are  examined  in  this  Appendix  with  E^  =  0  when 

the  larger  guide  is  filled  by  a  homogeneous  dielectric  (R^  =  =  0). 

Substituting  (52)  and  (53)  in  (56)  and  reauiring  that  E  .  =  0  relates  b  to  a 

xi  nm  nm 

as 

'  0  . 

b  =  -  a  — —  (99.) 

nm  nm  7 .  0  .  > 

1  xi 

B  is  related  similarly  to  A  by 
nm  nm 


B  =  A  ■ 

nm  nm  7  0 

P  xp 


After  introducing  the  notation 


(100) 


xil 

10  a10  !<%  7 

o  ll 


2  2 
P  .  +  0  . 

xi  yi 

anm  i^  7.0. 

01  xi 

2  2 

0+0 

»  xp  yp 

nm  i<^a  7  0 

0  p  xp 


(101) 


(102) 


(103) 


The  field  components  E  .,  H  .  and  H  .  may  be  computed  from  (52)  to  (55), 

y  j  ^  j  y  <3 

bl),  (58),  (59),  (99)  and  (100)  as 


E  .  =  i^p  •*—  6 . 

yj  o  oz  j 


H  .  =  (k.2  +-4)6. 

xj  a  cX2  J 


(104) 


(105) 


-’St 


H  .  = 

yj 


(106) 


where 


*y  2  _  'y  2 

6i  =  d10  (e  ^  '  Re  ^  )  sin  ^xilx' 


(107) 


V  "V 

+  )  ,  d  sin  3  .x1  cos  3  .y1  e 

nm  xi  yi 


D  sin  p  x  cos  0  ye 
ran  xp  yp 


(108) 


The  scalar  function  6 .  can  he  identified  as  the  x-component  of  a  magnetic 

J 

Hertz  vector  H.  =  ir  6.  ,  which  satisfies  the  wave  equation 
J  x  J 


sPn.  +  k  2n  .=  o 

-J  J  -  J 


(109) 


Tlie  electric  and  magnetic  fields  are  related  to  H,  as 

-J 


E .  =  i<up  V  x  n . 

-J  O  -j 


(110) 


h.  =  k  n.  +  v  (v  .n.) 

— j  j  -j  —  j 


(in) 


Continuity  of  the  tangential  electric  and  magnetic  field  components 
across  the  aperture  requires  that 


^  Si 


2-  6 
oz  p 


(112) 
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(113) 


(k£  +  )  6.  =  (k2  +  ~  )  6 

i  -n  2  i  v  -n  n.2/  p 

ox  ox 


dy  dx  ~  dy  dx 


(114) 


The  last  two  conditions  can  be  satisfied  for  a  y-dependent  6  .  only  if 

J 

k.  =  k  ,  which  is  the  trivial  case  of  no  dielectric  discontinuities.  However, 

i  p'  ’ 

there  are  no  contradictions  if  6^  docs  not  deperd  on  y  within  the  aperture. 

In  this  case  H  =  0,  H  does  not  depend  on  y,  and  the  continuity  of  the  tan¬ 
gential  field  components  requires  only  (112)  and  (113)  to  be  satisfied. 

In  a  cample tely  open  small  guide  (y_^  =  2e,  x^  =2 £  )  the  function  6^, 
which  does  not  depend  on  y,  can  be  represented  by  the  m  =  0  terms  of  the  double 
summation  (107).  Hence  E  .  (and  also  E  )  will  exhibit  no  y  variations  over  the 

yi  yp 

aperture .  Conversely  E  . ,  which  depends  on  y,  leads  to  a  y-dependent  6 .  and  to 

yi  i 

a  y-dependent  H^,  which  violates  the  continuity  of  H^.  A  y-dependent  E  and 

a  y-independent  H  .  lead  to  E  .  £  0.  (This  can  be  seen  from  (89),  (92)  and  (93) 

of  Aunendix  C  in  the  case  where  H  .  is  the  same  as  in  the  principal  wave  guide 

xi  0 

mode ) . 
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Figure  1  .  Waveguide  Geometry 
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of  the  Slot  Radiating  Into  a  Wide  Waveguide 


Figure  4.  Waveguide  Conductance  xq  =  yQ  =  10. 5x 
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Figure  10.  Variational  Approximation 

to  the  Voltage  Distribution  Along  the  Waveguide 
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(Figure  12.  Profile  of  Electron  Density  and  Collision  Frequency  for  o  Velocity  of  17.5  K 

FPS  ond  Altitude  of  200  K  FT 
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Figure  15.  Waveguide  Susceptance  Comparison  of  Solutions  for  E 


Figure  16.  Conductance  of  the  Waveguide  Slot 
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